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ABSTRACT 

Recent hydrodynamic simulations of galaxy formation reveal streams of cold (T < 



K) gas flowing into the centers of dark matter halos as massive as 10^^ ^'^■^Mq at red- 
shifts z ~ 1-3. In this paper we show that if >20% of the gravitational binding energy 
of the gas is radiated away, then the simulated cold flows are spatially extended Lya 
sources with luminosities, Lya line widths, and number densities that are comparable 
to those of observed Lya blobs. Furthermore, the filamentary structure of the cold 
flows can explain the wide range of observed Lya blob morphologies. Since the most 
massive halos form in dense environments, the association of Lya blobs with over- 
dense regions arise naturally. We argue that Lya blobs - even those which are clearly 
associated with starburst galaxies or quasars - provide direct observational support 
for the cold accretion mode of galaxies. We discuss various testable predictions of this 
association. 

Key words: cosmology: theory-galaxies: haloes-galaxies: formation-galaxies: cooling 
flows-galaxies: intergalactic medium-line: formation 



1 INTRODUCTION 

The physical origin of spatially extended Lya sources, also 
known as Lya blobs (LABs) , is still enatic. LABs have been 
associated with cooling radiation, in which gas that collapses 
inside a host dark matter halo releases a significant frac- 
tion of its gravitational binding energy in Lya line emission 
jHaiman et al.ll2000l: iFardal et allbOOll : iBirnboim fc Dekell 
l2003l : iDiikstra et all liood also see Katz & Gunn 1991). 
Other mechanism that have been invoked to explain the 
origin of LABs include: photoionization of cold (T ~ 10^ 
K), dense, spatially e xtended gas by o bscured quasars 
llHaiman fc Reed 1200 ll: iGeach et all 120091 '). population III 
stars ( Jimenez &: HaimanI bOOd). or spatially ex tended in- 
verse Compton X-ray emission (jScharf et al.ll2003l ) ; the com- 
pression of ambient g as by superwinds into den s e, cold Lya 
emitt ing shells (e.g. iTaniguchi fc Shioval |2000| : I Mori et ahl 
I2OO4I') ; or st ar formation that is triggered by relativistic jets 
l|Reesl ll989'l: or some combination of photoionization and 
cooUnglPurlanct to ct a l. 2005). 

Since their discovery by Steidel et al. (2000, also see 
Keel et al . 1999), several tens of new LABs have been 
found fe.g.lMatsuda et al.ll20oi : iDev et alllioosl : ISaito et all 



wavelength studies of Lya blobs have revealed a clear as- 
sociation of the brighter blobs with regular Lyman Break 
galaxies (LBGs, e.g. Matsuda et al. 2004), sub- millimeter 
(sub-mm) and infrared (IR) sources which imply star for- 



mation rates of 10^ Mq yr ^ ( Chapman et al] 1200 ll : 



I2OO6I : ISmith et all |2009^ . These are typically associated 



iGeach et al]|2005l. l2007l: iMatsuda et al.ll2007l). or^it h un- 
obscured ^Bunker et al.' I2OO3I; IWeidinser et al' '200|}_ and 
obscu red quasars (Basu-Zych & Scharf 2004; Gcac h et al.l 
I2OO7I : [Smith et all 120091 '). However, in other blobs this as- 
sociation has been ruled out, which has led to the conclu- 
sion that cooling radiation by cold accreting gas may have 
been observed (Nilsson ot al. 2006; Matsuda ot al. 2006 
Saitq_etal. 2008; Smith & Jarvis 2007; Smith ct al. 2008a 

We propose that the sources that are associated with 
the Lya blobs may have little to do in powering the spa- 
tially extended Lya emission. This proposal is physically 
motivated by recent hydrodynamical simulations of galaxy 
formation, which show that baryons assemble into galax- 
ies through a two-phase medium which contains filamen- 
tary streams of cold (T ^ lO'^ K) gas embedded within a 
hot gaseo us halo ( e.g 'Keres et al. II2OO5I: 

i 



Dekel fc BirnboimI 
Dck el et al.ll2009l : lKeres et al.ll2009l : 



with massive halos that reside in dense parts of t he Uni- 
verse (jSteidel et al.ll2000l : lMatsuda et al.ll20o3 . l2006l ). Muhi- 
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2006; Ocv irk et al.ll2 008: 

Agertz et al.ll2009l ) . The cold flow s contain ^ 5-50% of the 
total gas content (|Birnboim et al . 2007; K eres et al.ll2009h 
in halos with mas ses in the range Mhaio ~ lO^'^-lO^^'^M© 
(|Dekel et al.ll2009l ). As we will show in this paper, these 
cold streams are probable sources of spatially extended Lya 
emission. Under very reasonable assumptions we find that 
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the cold accretion model 'predicts' the existence of spatially 
extended Lya sources with properties reminiscent of the ob- 
served LABs around massive halos (Afhaio ^lO^^M©). 

The outline of this paper is as follows: we describe our 
model in §(2] In §[3] we present our results, before discussing 
our work in § [3] and presenting our conclusions in § [S] The 
cosmological parameter values used throughout our discus- 
sion are (Q.„^,nA,iib , h, as) = (0.27,0.73,0.046,0.70,0.82) 
l|Komatsu et al.ll2009t ). and we denote t he primordial heliu m 
abundance by mass as Ykc = 0.24 (e.g llzotov et al.lll997l ). 



2 THE MODEL 

Our model is based on the scenario for cold accretion tha t 
has emerged from recent simulations l|Keres et al.l |2009| ). 
'Hot' gas that is in hydrostatic equilibrium with the dark 
matter potential well at the virial temperature of the dark 
matter halo is in pressure eq uilibr ium with the cold gas (also 
see iFall fc Reedil985l : lReeslll989l) . which makes up a frac- 
tion /cold of the total gas mass of the halo, Mgas. We further 
assume that the gas mass fraction, /gas, is half of the uni- 
versal value. Hence, the total gas mass Mgas = /gasAfhaio = 
0.5Q,b/^MMhs.io = 0.08Mhaio. This accounts for a substan- 
tial fraction of baryons that may be locked up in (dwarf) 
galaxies that reside within the halo of interest. The as- 
sumed gas mass fraction /gas = 0.08 is in good agreement 
with the valu e derived for groups of galaxies ai z < 1 (e.g. 
iGiodini et al.l .2009). As the cold streams navigate to the 
center of their host dark matter halo, they are progressively 
heated by the release of a fraction of their gravitational po- 
tential energy through weak shocks l|Rees fc Ostrikerlll977l : 
iHaiman fc Recs 2001). This heating of the cold gas is bal- 
anced by radiative cooling , mostly through the Lya emis- 
sion line. Next, we describe the model in more detail. 

• We use the gas density profile fo r hot gas that is in 
hydrostatic equilibrium with an NFW l|Navarro et al.|[l997l ) 
dark mat ter potent i al we ll with a concentration parameter 
c = 3.8 (|Gao et a l.' 2008) at the virial temperature of the 
halo as derived by Makino et al (1998, their Eqs. 8 and 11). 
This profile is very similar to an isothermal /3-model, which - 
according to the simulations - p rovides an accurat e descrip- 
tion of the hot gas component (|Keres et al.ll2009h . In 14.61 
we show that our final results are robust against variations 
in the assumed density profile. 

• The hot gas reaches the virial temperature, rh(r) = 
Tvir = 1.9xlO'^K(Mhaio/lO^^M0)2''3[(l-Fa)/4], and co-exists 
in pressure equilibrium with the cold gas. Pressure equilib- 
rium between the hot ('h') and cold ('c') gas implies that 
nc{r)Tc{r) = nh(r)Tk(r), in which n(r) and T{r) denote the 
number density of particles and the gas temperature, respec- 
tively. We obtain Tc{r) under the assumption that cooling 
balances the heating of the cold fiows that occurs as they 
navigate into the center of the dark matter halo. 

We parametrize the gravitational heating rate by assum- 
ing that a fraction /grav (/hot) of the change in the gravi- 
tational potential energy of each gas element along its tra- 
jectory is converted into heat in the cold (hot) gas. The re- 
maining fraction (1 — /grav — /hot) is converted into additional 
kinetic energy of the gas element. Throughout the paper, we 
assume that the transfer of energy into the hot gas is neg- 
ligible, i.e. /hot ~ 0. This assumption appears reasonable 



given that the maj ority of the cold gas mass is in smo oth 
continuous streams (|Keres et al.ll2009l : iDekel et al.ll2009l ) , as 
opposed to di screte clouds w hich are likely to heat the hot 
gas (e.g. Dekel fc BirnboimI [200S . and references therein). 
Under this assumption , which is discussed in more detail in 
ij |4.6.2l /grav = 1 corresponds to infall at a constant velocity, 
while /grav = corresponds to free-fall. 

For the sake of simplicity, we adopt the conservative work- 
ing assumption of a constant non-zero /grav inside the virial 
radius (rvir) of galaxy halos and /grav = outside. The most 
recent Smoothed Partic le Hydrodynamical (SPH) simula- 
tions l|Keres et al.ll2009h indicate that the cold fiows prop- 
agate inward at an approximately constant speed, imply- 
ing /grav ~ 1; however. Adaptive Mesh Refinemen t (AMR) 
simulations (|Ocvirk et al.ll2008l : [Pekel et all 120091 ') indicate 
that the cold gas accelerates throughout its motion and that 
therefore /grav may be smaller 0. Future simulations might 
be able to refine our working assumption by resolving the 
precise dynamics and heating of the cold fiows. Hence, the 
heating rate per particle is given by 

\ J- GM{<r)iimp 
H{r) = /grav X \,2 '^ v{r)+H^{r), (1) 

where ^ is the mean molecular weight per particle in the cold 
fiow (in units of nip), and v{r) denotes the infall velocity 
which is given 

i;'(r) =2i;'(rvir) + 2(l-/grav-/hot) f ds ^^^^^ '\ (2) 

and we assume that f (rvir) ~ -\/2ucirc (with our results not 
being sensitive to this choice). The term H^(r) denotes the 
heating rate due to absorption of ionizing radiation. At the 
typical densities in the cold fiows (ric J>0.1 cm~'^) the gas 
is self-shielded from external ionizing radiation, and hence 
H^[r) = 0. It is possible however, that galaxies embed- 
ded within the cold fiow may photoionize some surround- 
ing re gion which would locally boost H{r) (|Furlanetto et all 
l2005h . In any case, ignoring this extra photoheating term 
only makes our predicted Lya luminosities from the cold 
fiows smaller, and therefore makes our results more conser- 
vative. As was mentioned above, we assume that /hot = 0. 

We obtain an equilibrium temperature Tdr) at radius r 
by equating H{r) to the cooling rate per particle in the cold 
fiow, which is given by, 

A(r, T, XHi) = — [nenHiCcooi(Tc) + nenmi-Rcooi(T'c)] . (3) 

Here, the first term in the square brackets denotes the 
cooling rate (in erg s~^ cm~'^) due to coUisional excita- 
tion of H atoms by electrons. The second term denotes 
the cooling rate due to recombination events of free elec- 
trons and protons (other cooling processes are negligibly 
small for the typical gas temperature in the cold fiow). The 
rate coefficients Ccooi(rc) and -Rcooi(7c) were taken from 
Hui & Gnedin (1997, hereafter HG97). Finally, the ioniza- 
tion state of the gas determines its cooling rate. Under the 

^ Our discussion focuses on the extended Lyo emission, and not 
on the compact core of the galaxy where the cold streams are 
finally brought to rest. Note that dust may suppress the Lya lu- 
minosity from the core, but is less likely to affect the cold streams 
which carry motal-poor intergalactic gas. 
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assumption that the gas is self-shielding (which is justi- 
fied later) we obtain a one-to-one relation between T and 
XHi through xm = arec,B (rc)/(arcc,B(rc) -I- Cion(rc)). Here, 
C'ion(Tc) denotes the coUisional ionization rate coefficient, 
and arcc,B(Tc) denotes the case-B recombination coefficient 
(as the cold neutral gas is optically thick in all Lyman-series 
lines and case-B applies; the related coefficients were taken 
from HG97). At the temperatures of interest, helium is neu- 
tral inside the cold gas and free electrons are only supplied 
by hydrogen, i.e. nm{r) = XHi(r)nc(r)/[l + (yHc/4)] and 
nHii(r) = ne(r) = [1 - a;Hi(r)]nc(r)/[l + (Yho/^)]. 

In practice we assume that Tc = 10* K and compute ric as- 
suming pressure equilibrium. Since temperature determines 
the ionization state of the gas in self-shielded gas, we ob- 
tain a cooling rate which we compare to the heating rate. 
If A(r, T, a;Hi) > H{r) {A{r,T,xm) < H(r)), then Tc is low- 
ered (increased) until the cooling and heating rates are equal 
to within 1%. 

• Once the temperature, density and ionization state of 
the gas have been determined, we compute the Lya emis- 
sivity (in erg s~^ cm~^) as a function of radius, 

<^hya{r) = nenmCLya{Tc) + 0.68huanennuarcc,B{Tc)- (4) 

Here, the first term denotes the luminosity density in Lya 
photons (in erg s~^ cm~'^) following coUisional excitation of 
H atoms by free electrons. The coUisional excitation coeffi- 
cient is gi ven by Clyq = 3.7 x 10~" exp{-hiya/kT)/T^^^ erg 
s~^ cm^ (|Osterbrockll 19891 . p 55). The second term denotes 
the luminosity density in Lya photons following case-B re- 
combination, and that ~ 68% of all case-B recombination 
events result (|Osterbrocklll989l ) in a Lya photon of energy 
hva = 10.2 eV. In practice, the recombination term can be 
safely ignored. 

• Lastly, we obtain the total Lya luminosity by integrat- 
ing over volume, namely 

iLya = / dr 47rr^eLya(r-)/coid^^ — . (5) 

Jo -'h 

Emission from r > rvir is expected to have a surface bright- 
ness that is well below the detection threshold of existing 
observations, because beyon d ryir no rarefied hot gas ex- 
ists t o confine the cold flow (iKeres et al.ll2009l : iDekel et al.l 
I2OO9I ). and so the cold gas density declines considerably. 
The gas at r > rvir therefore provides a negligible con- 
tribution to the total Lya luminosity. The fraction /com 
(which we as sume to be indepen dent of radius) is taken from 
Figure 3 of H eres et al] l|2009l) . which we parametrize as 
/cold = O.25(Mhaio/lOi2M0)-°-55 for Mhaio ^ 4 X 10"Mq, 
and /cold ~ otherwise. Furthermore, %, denotes the frac- 
tion of Lya that makes it to the observer. The factor 
/coid7c(r)/rh in the integrand denotes the fraction of the 
volume that is occupied by the cold gas in the range of radii 
within r ± dr/2. 



3 RESULTS 

Throughout the paper we assume that /grav ~ 0.3, which 
we regard as a conservative choice (see § [2] for a more de- 
tailed discussion). Furthermore, we assume that Ta = 0.5. 
That is, ~ 50% of the emitted Lya photons is observed. The 
most likely source of opacity is provided by residual HI in 
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Figure 1. Top panel: the solid and dotted lines show the radial 
dependence of the temperature and electron fraction by number 
of the cold gas phase for a dark matter halo of mass Mjjaio = 
IO^^Mq, and under the assumption that a fraction /grav = 0.3 of 
the gravitational binding energy of the gas is converted into heat. 
Central panel: the solid lines show the number densities of gas 
particles in the cold and hot phase and and the dotted line shows 
the infall velocity of the gas expressed in units of fcirc- Bottom 
panel: the fraction of Lya luminosity that is emitted interior to a 
radius r (normalized by the total, La^tot « 1.3 X 10*2 erg s"!). 

the intergalactic medium (IGM), which mostly affects the 
flux of photons that was emitted b lueward of the Lya res- 
onance (e.g Fig 3 of lMadau|[l995l ). At 2 = 2, z = 3 and 
z — 4, the mean transmitted fraction of photons emitted 
blueward of the Lya resonan ce is (%,) ^ 0.8, (Tg) ~ 0.7, an d 
{Ta) ~ 0.4, respectively (e.g. lFaucher-Giguere et al.|[2008bt ). 
However, these fractions apply to an average that was taken 
over a scales of tens of comoving Mpc (cMpc) along the line 
of sight. Since the massive halos of interest reside in dense 
parts of the Univ erse, we expect the local to be more opaque 
than average (e.g. lDiikstra et al.ll2007l : lDiikstrall2009l ). which 
motivates our choice Ta = 0.5 at z ^ 3. The gas in the cold 
flows consists mainly of compressed neutral intergalactic gas 
and dust opacity is likely negligible (see § 14.21 for additional 
discussion on dust). 

3.1 The Halo Mass-Lya Luminosity Correlation 

Consider a halo of mass A/haio ~ IO^^Mq, for which iicirc = 
220 km s"^ Tvir = 1.9 x 10^ K, and rvir = 82 kpc (see, e.g. 
Barkana & Loeb 2001, their Eqs. 24-26). The solid line in 
the top panel of Figure [T] shows the radial dependence of the 
gas temperature in the cold flow. The dotted line shows the 
electron fraction by number, i.e. Xe — xmi = nmi/rip, where 
Up denotes the number density of protons. The temperature 
is found to be in the narrow range of ~ 1-1.1 X 10* K which 
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Figure 2. Same as Fig [T] but for a dark matter halo of mass 
Mhalo = 1013Mq. 



results in a very low level of ionization {xe ~ 10~^). This 
is caused by the fact that atomic line cooling is very effi- 
cient at the densities encountered in the cold flow (see cen- 
tral panel), and so the gas cools down to the temperature 
floor below which collisional line excitation is ineffective. 
The temperature and ionization fraction increase outward 
as the density, and therefore the cooling rate decreases with 
radius. Figure [2] shows the same quantities for a halo of mass 

The solid lines in the central panels shows the number 
densities of gas particles in both the cold and hot phase. 
The dotted line shows the infall velocity in units of Vch-c- 
(Note that the infall velocity enters the heating rate through 
Eq. [T]) The cold gas is denser by a factor of Tvir/Tc- This re- 
sults in number densities in the cold phase that exceed 1 — 10 
cm~"^ at r ;^0.5rvir. These high densities justify our assump- 
tion of self-shielding: the recombination rate in the cold flow 
is arocWc ~ 10~^^ — 10~^^ s~^, while the photoionization 
rate by the UV ba ckground in the IGM at z = 3 is Fuv ~ 
6 X 10"^^ s"^ fe.g iFaucher-Giguere et~al]|2008al '). Even the 
gas on the very edge of the cold flow, that is optically thin to 
ionizing background, contains an ionized fraction of hydro- 
gen (by number) of only xui = 1 + ^ — |(ra^ -l-4a)^'^^, where 
a = Fuv/wcQrcc- Therefore, xui = 0.2 — 0.6 for ric = 1 — 10 
cm~'^. That is, even this gas is significantly neutral, and the 
gas becomes self-shielding less than a pc into the fiow. Note 
that boosting the local photoionization rate by a factor of 
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10^ does not change this conclusior0. 



Figure 3. The solid line shows the dependence of the total hya 
luminosity as a function of halo mass, Mhaioi when a fraction /grav 
of the gravitational binding energy of the cold gas is converted 
into heat. This model predicts Lya luminosities in the range 10*^— 
10*^ erg s~^ (comparable to the range of values observed from 
the majority of blobs, see Matsuda et al 2004) for halo masses in 
the range Mhaio = 10^^ — 3 X 10^^ Mq. To explain the luminosity 
of the most luminous (Llyq = 10** erg s~^) blob with cold ac- 
cretion requires a halo of mass M^aio ~ 3 X 10 ^^ Mp), close to th e 
maximum halo mass that can host cold flows llDekel et al.ll2009l) . 
The blue dotted line shows the soft X-ray (free-free) luminosity of 
the hot gas. There are not nearly enough X-rays to significantly 
boost the hya emissivity of the cold flows (see text). 



The lower panel shows the total Lya luminosity that 
was emitted interior to radius r (normalized to the total, 
which is L„,tot ~ 2 x 10*^ erg s'^). We find that ~ 50% of aU 
Lya is emitted at r ^O.Srvir- Hence, our total calculated Lya 
luminosity is emitted over a spatially extended region. It is 
not possible to robustly predict the surface brightness profile 
as this is determined by the number of filaments in the cold 
flow, their geometry and their orientation relative to the 
observer. The surface brightness profile would furthermore 
be sensitive to any radial dependence of the 'gravitational 
heating efficiency '-parameter /grav 

Figure [3] shows the dependence of the total Lya lu- 
minosity, Li^ya, as a function of halo mass, Afhaio as the 
solid line. We find Lya luminosities in the range 10*^ — 
10'*'^ erg s~^, compa rable to the observed range in LABs 
l|Matsuda et al.|[2004| j for halo masses in the range Afhaio ~ 
1O^^-1O"M0. In order to obtain LLya = 10*** erg s'\ which 
is the lumino s ity of the brightest of blobs that was found by 
ISteidel et al] (|200G| ). we either need a halo with a mass of 
~ 3 X IO'^^Mq, or a somewhat lower halo mass with a higher 

/grav • 

The blue dotted line shows the total X-ray luminosity of 
the hot phase which is given by, 

Lx= dr 47rr'ne(r-)nHii(r)(l + lHc)eff(TviO, (6) 

^0 

where eff(Tvir) = 1.4 x IQ-^'^T^/'^ erg s"^ cm^ is the free- 



^ Because in our model the Lya blobs are associated with mas- 
sive halos, Mjjaio > 10^^ M0, the cold filaments are compressed 
by a factor of (Tvir/lO* K) ;^2OO(Afhalo/lO^^M0)2/3 by the hot 
virialized gas. Therefore the gas in our model is denser and better 



capable of self shielding than in the models of iFurlanetto et al.l 
li2005i) . For this reason, the photoionizing background is not dom- 
inating the heating - and resultant Lyo cooling — rate from the 
cold flows 
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free emissivity (see e.g. Eq 5.15b of iRvbicki fc LightmanI 
[197^. Note that the hot gas is fully ionized with nB_ii{r) = 
nh{r)/{2 + ^Yhc), and an electron density Ue = (1 + 
|iHc)nHii. 

The X-ray luminosity is comparable to the Lya lumi- 
nosity. One may think that the X-rays can penetrate the 
cold flows and boost the Lya luminosity. This is unlikely the 
case, as the frequency dependence of the free-free spectrum 
is given by L^jf oc exp(— /ipi^/fcTh), and therefore the hot 
gas emits photons as energetic as hpU ~ 0.5(Tvir/5 x 10^ K) 
keV. These soft X-ray photons easily penetrate deep into 
the cold flows where they photoionize either hydrogen or 
helium atoms, and create an energetic electron. In a neu- 
tral medium, ~ 40% of the energy of the electron goes into 
exciting HI, ~ 30% goes into coUisionally ionizing HI, and 
~ 10% goes into heating the gas (the rest excites and ion- 
izes helium, ShuU & Van Steenberg 1985) . Hence, even if the 
cold flows absorb all X-rays, line excitation by secondary 
electrons would only boost the Lya emissivity by at most 
a factor of <;1.4. Furthermore, a boost in the ionization 
(and thus electron) fraction naively boosts atomic line cool- 
ing rate - and hence Lya emissivity - by the same factor. 
However, even if X-rays (significantly) enhance the electron 
fraction in the cold flow, the cooling at a given temperature 
would only be more efficient so as to balance the heating 
rate. Since at most 10% of the X-ray energy is converted 
into heat, this heating rate is boosted by a factor <1.1. 
The above estimates assume that the cold flows absorb all 
the X-rays; this is very unlikely to occur in reality because 
the cold medium occupies a fraction ~ fcoidTc/Ti^ ~ lO"**- 
10"'^ of the volume of the halo, and the majority of X-rays 
freely escape out from the halo. The boost in Lya emis- 
sivity from the cold flows due to X-rays is therefore neg- 
ligible. The X-ray photons that escape would redshift to 
energies hpi/ ^0.125(rvir/5 x 10® K) keV, and would not 
be detectable with existing X-ray telescopes. Note however, 
that metals can boost the X-ray emissivity of the hot gas at 
energies 3> 1 keV. Deep X-ray observations may therefore 
set useful limits on the amou nt of hot gas in the dark matter 
halos hosting the blobs (as in lGeach~e t al. 2009, who put an 
upper limit on the average mass of halos hosting Lya blobs 
of Mhaio < 10^^ Mq), although these limits do depend on 
the assumed metallicity of the hot gas. 



3.2 The Lya Luminosity Function of 
Cold- Accretion Powered Blobs 

The previously obtained relation between galaxy mass and 
Lya luminosity allows us to compute a luminosity func- 
tion of accretion powered LABs from the mass function of 
dark matter halos. Observed LABs are known to be asso- 
ciated with a region that contain a significant overdensities 
of Lyman-break galaxies (L BGs) with a fract i onal excess 
density of 5lbg = 5.0 ± 1.2 (|Steidel et al.lll998l . |2003| ). It is 
essential to take into account this overdensity when predict- 
ing tlM_ex£ected_jiur^ of LABs brighter than La 
as l|Barkana fc Loebll2004l ). 

?ibiob(> La\5M{Vs)) = dm — — B{m,SM{Vs)), 

(7) 
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Figure 4. The solid line shows the number density of Lya blobs 
brighter than La, «blob(> La) cMpc~^, and taking into account 
that fractional overdensity of the region of the Universe in which 
the blobs reside is Sm{Vb) = 3o"jvf(l4). The dashed and dotted 
line refer to values of (5m(K>) = 2crjv<f (Vs) and S^iVs) = 5(Tm{Vs), 
bracketing the likely range of overdensities inferred from the ob- 
served overdensity of LBGs in the region by Steidel et al. (1998). 
Here, <jj^,j{Vs) denotes the variance of the density field averaged 
over spheres of volume Vs = 1.3 X 10^ cMpc^ (the survey volume). 
Our model with /grav = 0.3 provides a good fit to the data. The 
present uncertainty on Sm{Vs) allows models with /grav ^0.15 
to reproduce the observed number of Lya blobs. The black dot- 
ted line shows the number of blobs expected in an 'average' part 
(i.e. Sm{Vb) = 0) of the Universe, indicating how rare the most 
luminous blobs are (nblob(> 10'*^erg s"-'^) ^ 10~^ cMpc~^). 

where {dnsT / dm) dm denotes the Sheth-Tormen number 
density of dark matter halos (in cMpc"'^, Sheth & Tormen 
1999) in the mass range Mhaio = m ± dm/2. The factor 
B{m,5M{Vs)) denotes the boost in the number of halos due 
to the overall matter overdensity, (5m (Vs), within the survey 
volume Vs. 

In Figure 13] we plot nbiob(> L^; (5m (Vs)) as a func- 
tion of Lya luminosity for three values of the overdensity 
Sm{Vs) = 2(jm(Vs) {dashed line), 3aM{Vs) {solid line), and 
4(Tm(Vs) {red dotted line). Here, o'm{Vs) = 0.013 denotes 
the variance in the matter density field at 2 = 3 averaged 
over spheres of the survey comoving volume Vs = 1.3 x 10^ 
cMpc^. The overdensities 5m{Vs) were chosen so as to re- 
produce the observed boost and standard deviation in the 
LBG number density, under the assumption that LBGs are 
associated with dark matter halos of mass A/haio ~ 10^^ Mq 
(see Appendix IX)) . 

The data points were d erived by taking the ratio be- 
tween the number of blobs (|Matsuda et al.l 120041 ) brighter 
than La and the survey volume V = 1-3 x 10^ cMpc'^. The 
error bars on the data points were obtained from the rela- 
tion (T^ = 6'^(Tm(V) + o-p, where ap denotes the Poisson 
error and b^a%[{VB) denotes the cosmic variance for a bias 
fac tor 6 = b(Mhaio, z ), wh ich is in turn calculated follow- 
ing [Somerville^^et^l] (|2004l ). For the cumulative luminosity 
function, the errors were added in quadrature. 

Figure |4] shows that the model luminos ity function 
agrees well with the data (|Matsuda et al.|[2004l ). The uncer- 
tainty in (5m(Vs) allows for a large spread in the predicted 
"-biob(> La), especially at the bright end. This is because 
the most massive (Mhaio ~ W^^Mq) halos are highly bi- 
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ased tracers of the mass and their number density depends 
sensitively on Sm{Vs). The sensitivity of the cumulative lu- 
minosity function to 5m{Vs) implies that good fits to the 
data are possible for a range of gravitational heating efficien- 
cies, /grav In the allowed range between Sm{Vb) = 2(TA/(Ki) 
and AaniiVs), we obtain a good fit for /grav > 0.15. Our 
predicted number of Lya blobs is practically degenerate in 
the product x /grav, and this constraint can be recast 
as Xi/grav >0.075. Therefore, only a relatively low gravita- 
tional heating efficiency suffices to make the cold flows 'glow' 
in Lya at luminosities observed from the Lya blobs. 

Lastly, the black dotted line in Figure 3] shows the num- 
ber density of blobs expected in an 'average' part (i.e. 
SMiVs) = 0) of the Universe. This implies that the most 
luminous blobs that were originally discovered by Stei- 
del et al (2000), are significantly rarer in the field (where 
Tibiob(> lO^^erg s~^) ~ 10^^ cMpc~^). Our prediction is 
consistent with a recent upper limit of nbi ob(J/a > 10^* erg 
s"^) ^5 X 10"'' cMpc"^ (|Yang et al.ll2009l '). 



4 DISCUSSION 

4.1 Predicted Lya Blobs Spectral Line Shapes 
and Physical Sizes 

The t ypical velocity widths (FWHM) of the LABs span the 
range ijMatsuda et al.ll2006l ') of 500-1700 km s'^ In the cold- 
accretion model, infall occurs at velocities in the range 1.0- 
2.5«circ ^ 220-1500 km s~^, which naturally result in the 
observed velocity widths (taking into account projection ef- 
fects due to the orientation of the filaments, and the inter- 
galactic absorption that suppresses the observed fiux of the 
'bluest' Lya photons (e.g. iMadau.1995. )). Although the gas 
in the cold flows is mostly neutral, resonant scattering is 
not expected to broaden the Lya line beyond the observed 
velocity widths of the blobs, as we discuss next. 

The cold flows are highly elongated structures, and the 
majority of the Lya photons escape in a direction perpen- 
dicular to the axis of the flow. Therefore, at the densities of 
interest the Lya photons have to traverse an HI column den- 
sity iVni ~ 10^^ - 10^^(i?flow/0.5 kpc) cm before escaping 
from the flow. Here, Rnovr denotes the radius associated with 
the characteristic cross-sectional area of each streanjfj. Scat- 
tering through such large columns broadens the Lya line to 
a width of ~ 800(jVHi/10^^ cm-2)^''^(Tc/10* K)^/^ km s"^ 
(|Neufeldlll99Q V Velocity gradients that exist throughout the 
flow will reduce this broadening: when velocity gradients 
are present, the relevant length scale - and hence column 

^ The cross-sectional radius of a cold filament can be estimated 
as follows. The total solid angle that is covered by the cold 
gas at a distance r from the center of the halo is SlcoidC'") = 
4tt X (/cold;^^)- Cold accretion typically occurs via multiple cold 
streams, and the average solid angle that is covered by a single 
cold filament is r2fii(r) = X 47r X (/cold^^)i where Nfn is the 
total number of filaments. From this it follows that the angular 
size of a single cold flow in radians, when viewed from the center 
of the halo, is Bai = ^ X 4 X (/cold^^)i which translates to 
a physical radius of i?flow('') = '^fll''- Substituting numbers yields 
^flow = 0.8(Affii/5)~^/^(r/0.5rvir) kpc (where the mass depen- 
dence enters only through /cold )■ 



density - that set the amount of velocity broadening is the 
Sobolev length, defined as the length over which the veloc- 
it y changes by one thermal width Is = Vth/[dv{r, R)/dR] 
(iBonilha et al.lll979l . and see e.g. § 6.2 in Dijkstra & Loeb 
2008b). Here, R denotes the transverse direction in the 
cold fiow. If for example, dv/dR ~ Vdrc/R (which appears 
reasonable for some cold filaments based on inspection of 
Figs. 8 & 9 of Dekel et al 2009), then we get that the to- 
tal amount of velocity broadening is reduced by a factor of 
(t^circ/wth)'''^ ~ 2.6(A/haio/lOi2M0)i/^ We therefore con- 
clude that although scattering plays a non-negligible role in 
broadening the Lya line profile, it does not lead to Lya ve- 
locity widths that exceed those of the observed blobs. How- 
ever, the complexity of the associated radiative transfer ef- 
fects does not allow us to make precise predictions for the 
Lya line width and its relation to the halo mass (or corre- 
sponding Lya luminosity). 

Because of the low volume filling factor of the cold fiows, 
we do not expect the Lya photons to scatter off neutral hy- 
drogen (HI) in separate cold filaments once they have es- 
caped from their filament of origin. For this reason, scat- 
tering will lead to a very low level of polarization (or none 
at all) for these models, in difference from the models of 
iDiikstra fc Loebl (|2008l ) who predicted high levels of polar- 
ization for Lya emerging from neutral, spherically symmet- 
ric, collapsing gas clouds. In addition, we do not expect 
the systematic blues hift of the Lya line that was predicted 
IDiikstra et al.l (|2006l ) for spherically symmetric modelfl 

Some resonant scattering would occur for photons that 
enter the IGM blueward of the line resonance. The impact of 
the IGM would be frequency- dependent, and we expect the 
associated 'blurring' to be most prominent for the most ener- 
getic Lya photons. The resonantly scattered Lya radiation 
is expected to have a lo w level ( ^7%) of linear polarization 
IDiikstra fc Loebll2008l ). 

The observed LABs have a large spatial extent 
ISteidel et al.ll2000D of ~ 150 kpc. In our simplified treat- 
ment, Lya is emitted from a region that extends out to the 
virial radius of the host halo. 

This implies LABs with diameters of ~ 160-350 kpc for halo 
masses in the range of Mhaio ~ IO^'^-IO^'^Mq. This range is 
expected to exceed the observed blob sizes, since our es- 
timate was derived under the simplifying assumption of a 
gravitational heating efficiency /grav that is independent of 

Interestingly, lAdams et aP l|2009^ did recently detect a system- 
atic blueshift of the Lyo line relative to the observed 21-cm ab- 
sorption line. They simultaneously reproduced both observations 
with radiative transfer of Lyo through a spherically symmetric, 
neutral collapsing gas cloud with a gas mass of Mgas ~ 10^^ Mq. 
Unless the mass of the host dark matter halo is implausibly high, 
A-^halo ~ lOl'*iW0, this would require an unusually high fraction 
of the gas to be in the cold phase (with /cold close to unity) . This 
inference might imply a significant scatter in /cold for a given halo 
mass. Some scatter in /cold is indeed observed in simulations, al- 
though not at the level to explain cold fractions close to unity. 
The possible existence of large /cold values would only boost the 
predicted number of cooling powered blobs and strengthen our 
basic conclusions. 
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radius . Higher resolution AMR simulations bv lAgertz et alj 
(|2009l ) show that cold streams slow down progressively after 
they enter the inner halo (r ^0.3 — O.Srvir), and therefore 
suggest that /grav is likely to increase at smaller radii. Mod- 
els that account for this increase predict smaller blob sizes 
that are closer to the observed values. It is also possible that 
existing observations only detect Lya radiation from the in- 
ner regions of cold flows where the Lya surface brightness is 
high. Indeed, ~ 60 — 70% of our computed Lya luminosity 
emerges from r < O.Srvir (with this modest reduction factor 
having a negligible effect on Fig. [Jjl . In this case, deeper ob- 
servations would reveal filamentary extensions of the blobs 
to larger scales. 



4.2 The Role of Dust 

Neutral gas in the cold flows forces each Lya photon to scat- 
ter multiple times before escaping from the flow. As a result, 
Lya photons traverse on average a total distance through 
the cold flow that is ~ 130(iVHi/10^2 cm-^) (1^/10^ KW'' 
times larger than that for continuum photons (|Adamslll975h . 
One would therefore expect the Lya flux to be highly sus- 
ceptible to absorption by dust. In reality, dust extinction is 
probably not an important effect since the cold flows consist 
of compressed intergalactic gas that has not been processed 
through galaxies. We therefore expect these flows to con- 
tain significantly less dust than the interstellar medium of 
star forming galaxies. Indeed, it is likely that quenching of 
Lya fiux by dust is more important for superwind-generated 
emission (see § 14.41 for a discussion of alternative models for 
Lya blobs), since in this case the Lya flux is generated in 
cold neutral shells of gas that we re swept up from the in- 
terstellar medium of the galaxy. iMori et al.l ((2004) found 
the dust-opacity to Lya of their superwind generated shells 
to be negligible (r^ ~ 0.03[Nm/5 x 10^^ cm"^]). However, 
their estimate of Td does not take into account the fact 
that scattering can boost the total distance traversed by 
Lya photons through the shell - and therefore the opti- 
cal depth-by a factor of ~ 110(Tc/10* K)^/*^ to a value of 
Ta ~ 3.3{Nm/5 x 10^^ cm-2)(Tc/10* K)^^^. It is therefore 
a serious concern that the dust opacity to Lya photons in 
these models exceeds unity. 

On the other hand, the cold flows are expected to con- 
tain significantly less dust because they originate in the , and 
so we conclude that the Lya flux generated in cold flows is 
not subject to the same level of dust opacity. For a ty pical 
inter galactic gas metallicity Z ~ 10"^^© (e.g. [Schave et al.l 
and for typical HI column densities and velocity gra- 
dients in the cold flows (see § 14. ip , we expect a dust opacity 
that is well below unity. This low level of enrichment is not 
unreasonable, as galactic outflows would typically avoid the 
overdense inflowing fllaments. Indeed, preliminary simula- 
tions indicate that outflows generally enrich the large volume 
surrounding the voids, while leaving the cold flows metal 
poor (R. Joung, private communication). 



to be surrounded by cold flows. Simulations support this 
assumption: the massive halos which produce a detectable 
Lya luminosity typically reside in overdense regions of the 
Universe in which accretion of cold gas from filaments onto 
the more massive halos is a continuous process at redshifts 
z > 2. This is in sharp contrast to intense starburst and/or 
quasar activity, which occur with duty cycles that are much 
smaller than unity. Indeed, in our model the cold fiows are 
luminous in Lya irrespective of the activity in the central 
galaxy, which explains why only a fraction of Lya blobs are 
observed to be associated with intense starburst or quasar 
activities. 



4.4 Comparison to Other 'Blob Models' 

While the cold accretion model can successfully reproduce 
the abundance of LABs in the Universe as well as their phys- 
ical properties, it is useful to compare our model to the other 
models for LABs mentioned in § [T] The observed number 
density of blobs, combined with their preferred residence in 
dense regions of our Universe, requires an association with 
massive halos (Afhaio ~ 10^^-10^^ Mq). 

An association with bri ght quasars also natu rally places 
the blobs in massive halos (jHopkins et al.|[2007^ . However, 
the filamentary geometry of the simulated cold gas shield it 
from any central source of ionizing radiation, suppressing the 
photoionization rate by the q uasar. Furthe rmore, the short 
duty cycle of quasar activity (|Martinill2004h is problematic. 
Indeed, a large frac tion of blobs have no associated quasar 
nearby (sec Bunke r et all l2003l : [WcidinEcr et al. 2004, for 
examples of Lya blobs associated with radio quiet type I 
Active Galactic Nuclei (AGN)). 

The duty cycle argument, combined with the lack of 
the relevant associated sources, argues against photoioniza- 
tion by inverse Compton radiation, or star formation that is 
triggered by AGN jets. Formally, photoionization by regular 
star forming galaxi es (LBGs) is not rule d out for ~ 60% 
of Lya blobs (e.g. iMatsuda et al.l l2004h . However, it re- 
mains to be explained why the LBGs in these dense en- 
vironments would be so efficient at generating Lya radia- 
tion, while observationally it is known that the escape frac- 
tion of Lya radiation from LBGs is typically low, and that 
only 20 — 25% of g = 3 LBGs would even classify as a Lya 
emitting galax y (Shaplev et al.l [20031 ) . Indeed, for this rea- 
son Jimenez fc HaimanI (|2006l ') argued for photoionization 



4.3 The 'Duty Cycle' of the Cold Accretion Mode 

When computing the cumulative luminosity function of 
blobs through Eq. (2), we implicitly assumed a duty cy- 
cle of unity. That is, we assumed each massive galaxy halo 



by primordial galaxies which can emit more ionizing pho- 
tons per observed rest-frame UV continuum. This is an in- 
teresting suggestion, although it remains to be shown that 
the gas in highly overdense, evolved, regions of our Universe 
can remain pristine down to 2; = 3 at a level that is needed 
to explain the Lya blobs. 

In any case, all photoionization models require cold spa- 
tially extended gas to be present in these halos, because the 
gas at the virial temperature of these halos (Tvir ^2 x 10^ K) 
would recombine too slowly to reproduce the observed Lya 
luminosities. Indeed, iHaiman fc ReesI l|200lh obtain their 
Lya recombination radiation - following photoionization by 
the central quasar - from cold (T = lO** K) gas clouds em- 
bedded in a hot medium, very similar to the cold-accretion 
scenario that is seen in galaxy-formation simulations (the 
main difference being that the cold phase consists of sepa- 
rate clouds in the early models [such as in Fall & Rees 1985] 
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rather than long continuous streams). Because photoioniza- 
tion models require the existence of cold spatially extended 
gas, they may be viewed as a special case of the cold accre- 
tion model in which the heating rate H{r) is dominated by 
photoheating (although it remains to be shown that the gas 
does not self-shield). In principle the photoionization mod- 
els also work when the cold gas is supplied by outflowing 
material (see below). However, dust in outflow models is a 
bigger problem here than for the inflow models (see § 14. 2p . 

A fundamentally different model for blobs is the su- 
perwind model. The association of some blobs with sub- 
mm sources imply a connection with starburst galaxies 
(having star formation rates of ~ 10^~^Mq yr~^). The 
kinetic energy associated with the starburst-driven super- 
winds is sufficient to power the Lyg emiss ion in LABs 
l|Taniguchi fc Shioval [2OO0I : iGeach et al.ll2005h . However, it 
is not clear that this coupling of energy occurs naturally at 
the le vel that is required to explain the LABs. I Mori et al.l 
l|2004 ) used hydrodynamical simulations to show that su- 
perwinds driven by a starburst (M«,max ~ 200 Mp, yr~^) can 
produ ce an LAB that closely resembles one of ISteidel et al.l 
l|200d ) blobs in appearance, but with a luminosity that is 
~ 10/Ta times too small (note that Ta< 1). If we simply 
scale up the luminosity of an LAB in proportion to the to- 
tal star formation rate, then the observed Lya luminosity 
requires a star formation rate of ~ 2000/7^ M© yr~^. How- 
ever, this is close to the observed star formation rate derived 
from the dete cted sub-millimete r flux in only one of Steidel 
et aL's LABs lGeach et al.ll2005h . Also, the conversion from 
star formation rate to the blob's Lya luminosity assumed 
that the dust opacity for Lya photons is negligible. We 
showed in I4.2l that this may not be correct when radiative 
transfer effects are accounted for. Furthermore, not all blobs 



have starb urst activity assoc iated with the m fMatsu da et al 
20061 : iNUs son ct al. 2006; S mith fc Jarvisl ^2007: Sa ito et al 



20081 : ISmith et al.ll2009l ') . thus ruling out this model for these 



4.5 Testable Predictions of the Cold Accretion 
Model For Lya Blobs 

Simulations suggests that the cold flows enter the center of 
the halo as multiple fllamentary streams. This filamentary 
structure is expected to show up in the observations. In §[4T] 
we argued that some 'blurring' of the Lya image may be 
expected, but only for the most energetic Lya photons, as 
these are expected to resonantly scatter in the IGM. Existing 
images of LABs typically have an angular resolution of ~ 1 
arcsec, which corresponds to ~ 8 kpc at 2: = 3. Therefore, 
existing observations would not resolve the filaments yet, al- 
though the filamentary nature of the cold fiows may show up 
at larger radii (where different filaments are well-separated). 
Existing observations of blobs that are presently thought to 
be powered by cold accreti on do show some interesting ir- 
regula rities in their images (jNilsson et al ] |2006l : ISmith" et al.l 
l2008al 'l. but the quality of the data does not allow to classify 
the image as 'filamentary'. It is intriguing that th e ima ges of 
some of the LABs discovered by Matsu da et al.l (|2004i) are 
very irregular, and possibly filamentary (e.g. blob 6, 9, 11, 
and 12). If these blobs are associated with cold flows, then 
we expect deeper and higher resolution images of the blobs 
to more clearly reveal the filamentary structure. 



We argued in § 14.11 that we do not expect the system- 
atic blueshift of the Lya line, nor the high level of polar- 
ization, that were predicted by spherically symmetric mod- 
els of neutral collapsing gas cloud. This is mostly because 
the Lya photons scatter mostly in the filament where they 
were produced, after which they escape to the observer (fol- 
lowing some resonant scattering in the IGM). Because the 
fraction of cold accreting gas increases towards lower galaxy 
masses (and therefore lower blob luminosities in our model) , 
we may expect cold accretion to be less fllamentary at lower 
galaxy masses. This implies that the predictions of the mod- 
els that assumed spherically symmetric accretion - involving 
the blueshift and the high levels of linear polarization - may 
be increasingly relevant at lower halo masses and blob lumi- 
nosities. However, these trends are rather weak in the limited 
range of halo masses of interest here, Mhaio = 10^^-10^'^ M© 
(see e.g. Fig 5 and Fig 6 of Keres et al. 2009, but see Adams 
et al. 2009 and ^WJ}. 

In our model, the vast majority of the Lya radiation is 
generated from coUisional excitation of atomic hydrogen. At 
temperatures of Tc ~ 10* K, coUisional excitation of the Ly/3 
transition is ~ 10 times less effective. Therefore, we expect 
an Ha (A = 6536 A) ffux that is ~ lOT^ x 6536/1216 ~ 54T^ 
times lower. This is unfortunate because a combined mea- 
surement of the Lya and Ha fluxes could have constrained 
the importance of radiative transfer effects. In 14.11 we con- 
cluded that radiative transfer effects were not negligible, and 
that they may noticeably broaden the Lya line. To find this 
broadening requires the detection of the much fainter Ha 
line. It is not guaranteed that the Ha flux will actually be 
suppressed by a factor as large as that mentioned above. 
As we argued in § 14.61 the weak shocks that heat the cold 
flows on their way to the halo center, may occasionally heat 
the gas to a temperature at which hydrogen becomes signifl- 
cantly ionized. In this case, the Ha flux may be down by the 
value expected from recombinations (~ 9Ta). These shocks 
may also result in helium becomi ng singly ionized, leading 
to a detectable He 1640 emission (|Yang et al.ll2006l ). 

The observed LAB luminosity function (which we c om- 
piled from the data presented by Matsuda et al.l l|2004l )) is 
reproduced by models in which 7^/grav ~ 0.15 (for SmCVs) = 
3(TAf(Ki), see Fig. 1). Assuming that 7^/grav ~ 0.15 (inde- 
pendent of redshift), we can predict Lya luminosity func- 
tions for accretion powered LABs at 2 = 2 and 2 = 4 
(when averaged over large regions of our Universe, such 
that 5m{Vs) ~ 0). These luminosity functions are shown 
in Figure 5, which implies that a decline in the number den- 
sity of accretion powered blobs more luminous than a few 
xlO*^ erg s~^ is expected beyond 2 = 3. Such a decline ~ if 
present - could be detected with future observations. For ex- 
ample, the Hobby-Eberly Te lescope Dark Energy Experiment 
CHETDEX. lHiU et al.|[20o3 ') is expected to detect 0.8 million 
Lya emitting galaxies and therefore thousands of LABs at 
1.9 < z < 3.5. Such a large sample of blobs will allow for 
accurate measurements of their spatial clustering, and lu- 
minosity functions. The combination of these measurements 
should place tight constraints on blob formation models. 

4.6 Discussion of Model Uncertainties 

We have taken the total fraction of the gas that is in cold 
streams as a function of halo mass from the most recent 
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Figure 5. The predicted cumulative luminosity functions of cold- 
accretion powered LABs at 2 = 2 (red dotted line), 2 = 3 (solid 
line), and 2 = 4 (blue-dashed line) (in a region at mean cosmic 
density) for /grav7Q= 0.15 (the value preferred by fitting to the 
data of Matsuda et al 2004; see text). A drop in the number 
density of sources more luminous than a few xlO'*'' erg s~^ is 
expected beyond 2 = 3. This drop may be enhanced by the that 
becomes increasingly opaque with increasing redshift. 

simulations of cold streams in massive halos (see § [l]) . Other 
model assumptions that were described in §[2] were based on 
these same simulations. However, the cold flow properties 
in our model - in particular the gas density and tempera- 
ture - can be quite different from what is seen in the sim- 
ulations. These differences are mostly due to the simplified 
nature of our model which assumes a two-temperature bary- 
onic component, while in reality there is a wider spread in 
temperatures in the cold streams (a more detailed discussion 
follows below). We have not attempted to bring our model in 
closer agreement to the simulations because it would require 
more detailed modeling of increasingly uncertain physical 
processes. Most importantly, we will show below that these 
differences are not important as far the Lya emissivity of the 
cold gas is concerned. In the following we address the dif- 
ferences in cold flow properties between our model and the 
simulations, and discuss other model uncertainties in more 
detail. 

4.6.1 Our Model versus Simulations: Cold Flow 
Temperature and Density 

Our models predict gas densities that are signific antly higher 
than those observed in AMR simulations (e.g. iDekel et al.l 
I2OO9I '). For a Mhaio = 10^^ Mq halo, we find ~ 1 cm"^ 
at r = O.Srvir, which is about a factor of ~ 30 larger than 
the mean cold fiow density that is found in the simulations. 
Similarly, the hot gas in our model denser by about a factor 
of ~ 3 — 4 (A. Dekel, private communication). There are 
several reasons for this discrepancy: (i) the simulated flows 
are at a somewhat lower redshift (z — 2.5) than we consider 
(z ~ 3), which reduces the density by a factor of (4/3.5)^ ~ 
1.5; (a) the cold fraction of gas is slightly higher in AMR 
simulations (/coM.amr = 0.4-0.5 at Mhaio = 10^^ M©, see 
Birnboim et al. 2007) than in the SPH simulations, which 
reduces the number density of particles in the hot phase. 
This amounts to an additional factor of (1 — /coid,SPH)/(l — 



/cold, AMR.) ^ 1.5. When combined, this explains a factor 
of (1.5)^ = 2.25. The remaining difference in the hot gas 
density is at the factor ^2 level, which may be due to the 
details of gas density profile that was used (also see i; l4.6.2|) . 

The fact that the density discrepancy is stronger for 
the cold gas suggests that pressure equilibrium is not ex- 
actly satisfied in the simulation. Indeed, in AMR simula- 
tions Pcoid ~ 0.5Phot at r <;0.5rvir (A. Dekel, private com- 
munication) which implies that at a given temperature, the 
number density in the simulated cold gas is reduced by an 
extra factor of ~ 2. Furthermore, the gas temperatures in 
the simulated cold gas are somewhat higher (a more detailed 
discussion of this follows below) which can lower the density 
by an extra factor of ~ 2. When collecting all factors of 2 
and 1.5, an order of magnitude difference in gas density can 
be accounted for. The remaining difference in the cold gas 
density is at the factor ~ 3 level, which may again be due 
to the details of gas density profile that was used (also see 
§|4X2]). 

We emphasis that because of our requirement that the 
cooling rate balances the local heating rate, H(r), our pre- 
dicted Lya luminosity is not affected by the exact density in 
the cold gas at all. This is because the heating rate (Eq[T} is 
independent of gas density and temperature. Furthermore, 
radiative cooling is dominated by the emission of Lya pho- 
tons at the typical cold gas temperatures. We therefore in- 
directly regulate the Lya emission rate per hydrogen atom 
by the local heating term. When the cold flow gas density 
is lowered/raised, a flxed cooling rate is maintained by rais- 
ing/lowering the cold gas temperature. For a fixed cold gas 
mass, we therefore predict a fixed Lya luminosity that de- 
pends almost entirely on the heating rate. Formally, there is 
a (very) weak temperature -and hence density- dependence 
of what fraction of the total cooling (and therefore heat- 
ing) rate emerges as Lya line photons. For example, at the 
temperatures that we encounter in our model ~ 63 — 67% 
of the all cooling radiation is in Lya line photons, while at 
T=2 X 10* K this fraction reduces to ~ 50%. The resulting 
very weak density dependence of the predicted Lya lumi- 
nosity has no impact on our results. 

In our model, the halo gas consists of 'hot' and 'cold' 
components. The hot gas is at virial temperature of the halo, 

Tvir = 1.9 X lO'^K(Mhalo/lO^^M0)2/3|(;^ _^ ^y^^ ^Jjjjg ^J^g 

cold flows have a typical temperature of Tc ~ 0.95-1.2 x 10* 
K due to efficient cooling (as the atomic cooling rate of the 
gas is high above lO^K and drops sharply below this value). 
Simulations show a wider spread in temperatures of the cold 
gas, namely Tc ~ 10*-10^ K. This apparent discrepancy is 
not serious, as the gas at 10^ K (in the SPH simulations) 
is mostly concentrated near the outermost regions of the 
halo where cooling is less efficient. The vast majority of the 
gas in the inner regions is in the range T ~ 1 — 2 x 10* 
K. Furthermore, the SPH simulations do not take into ac- 
count self-shielding of the gas (D. Keres, private communi- 
cation), which implies that photoheating is boosted artifi- 
cially. More importantly, ignoring self-shielding artificially 
suppresses the efficiency with which the gas can cool (e.g. 
see Fig 2 of Katz et al. 1996). This at least partially ex- 
plains the somewhat higher temperatures observed in these 
simulations. 

Lastly, the processes that are responsible for heating 
the cold flows may occasionally heat the gas to higher tem- 
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peratures, T ~ 2 x 10*-10^ K, after which it would rapidly 
cool to lower temperatures. The warmes t regions could b e 
bright in emission of the He 1640 A line ijYang et al.ll2006h . 
but at levels that are at most ~ 1-2 orders of magnitude 
lower than the total flux in the HI Lya line (which is barely 
affected by these rare warmer regions inside the cold flow). 
Hydrogen is mostly ionized at T ~ 2 x 10* K, which results 
in recombination radiation. This would boost the expected 
Ha flux emitted by the cold flows relative to the Ha-flux 
expected from purely coUisionally excited HI (see § I4.5|l . 

An important caveat to the simulations of cold flows is 
that they do not incorporate superwinds that are generated 
by star-forming regions. Superwinds can drive large-scale 
outfl ows that act as a source of Lya emission (|Mori et al.l 
l2004l ). These winds expand most efficiently into the lower 
density regions and will probably not affect the cold flows 
because of the small solid angle occupied by the cold flla- 
ments. Indeed, preliminary analysis of simulations that do 
take into account winds support this claim (D. Keres, private 
communication). Therefore, winds are expected to provide 
an additional source of Lya radiation, but would not dis- 
rupt the formation of Lya photons by gravitational heating 
of cold streams. Below we show that this additional 'wind- 
generated' Lya flux is likely comparable or less than the flux 
from the cold flows. 



4-6.2 Energy Dissipation in the Hot Gas and Dynamical 
Instabilities in the Cold Flows 

We assumed that no energy was dissipated in the hot gas, 
i-e- /hot = in EqO This assumption was motivated by the 
simulations in which the majority of the cold gas mass is in 
diffuse continuous streams. This implies that the majority of 
cold gas is surrounded by other cold gas, and that therefore 
little interaction occurs between the hot and the majority 
of the cold gas. However, shearing motions at the interface 
of the cold and hot gas may trigger the formation of insta- 
bilities that are not captured by existing simulations of cold 
flows. In particular, Kelvin-Helmholtz instabilities that de- 
velop at the cold-hot gas interface could completely disrupt 
the cold strearr[3, but this requires the relative velocity of the 
cold and hot gas to be subsonic i n both frames (e.g. iMiled 
I1957I : iBassett fc WoodwardI 1 19951 . and references therein) . 
Because the sound speed in the cold gas is Cs,coid ~ 13 km 
s~^, this requirement is clearly not met. However, the forma- 
tion of instabilities in cold streams is clearly an issue that 
needs to be addressed in future work. Instead of disrupt- 
ing the cold streams, these instabilities can introduce non- 
negligible amounts of energy dissipation in the hot ga^. 

We investigate the effect of energy dissipation in the hot 
gas in Figure |6] Here, we plot Lya luminosity as a function 
of halo mass for a model a model in which /hot = 0.7, i.e. 
~ 70% of the gravitational work done on an H-atom goes 
into heating the hot gas, instead of going into bulk infall 
motion of the atom (Eq[2)). Furthermore, for completeness 

® http://mngrid.ucsd.edu/ akritsuk/renzo/kelvin/kelvin.html 
^ Interestingly, these same instabilities would heat the cold gas 
through weak shocks, thus providing another physical mechanism 
by which gravitational binding energy is converted into heat in 
the cold gas. 



"AMR model" (text) 




Figure 6. The same as Fig [3] Here, the red dashed line repre- 
sent a model in which a significant fraction of the binding energy 
is dissipated in the hot gas, /hot = 0.7. Furthermore, we mod- 
ified densities in the hot and cold gas to more closely resemble 
those seen in AMR simulations (see text). The overall predicted 
Lya luminosities are reduced by a factor of ~ 2.. This can be 
compensated for by increasing /grav to /grav — 0.5. 



we adjust the cold flow properties to more closely resem- 
ble those encountered in AMR simulations. In this 'AMR 
model' we(i) assume different hot gas density and tempera- 
ture profiles as described in § 2.1 of Dekel & Birnboim (2008, 
their model with Ug =0), which appear to provide a better 
description of the hot gas component in the AMR simu- 
lation, (a) reduce the overall gas density by reducing the 
total gas mass fraction to /gas = 0.05, (Hi) boost the overall 
fraction of the gas in the cold phase by a factor of 2, (iv) 
assume that Pcoid = 0.5Phot. Modifications (ii-iv) reduce the 
gas density in both the cold and hot phase, and slightly in- 
crease the gas temperature in the cold gas (but not to the 
values that are encountered in the simulations, see I4.6.ip . 
The overall reduction in the amount of gas by a factor of 
~ 0.08/0.05 =~ 1.6 is compensated for by the boost in the 
cold fraction. However, the reduced infall velocity reduces 
the local heating rate and thus the overall Lya luminosity 
for a given halo. In order to restore the original halo mass- 
Lya luminosity relation, we require that /grav = 0.5. Note 
that it is possible that /grav + /hot > 1, and it corresponds 
to a model in which the cold gas decellerates as it descends 
down the gravitational potential well. In order for solutions 
to exist for v{r) at all radii, we require that /grav+/hot <1.35. 
Therefore a significant of the gravitational binding energy 
may be dissipated in the hot gas without invalidating our 
model, and for the example discussed above /hot <0.85 is 
tolerated. 



4.7 Comparison to Previous Work 

The high density gas in the cold flows can cool very ef- 
flciently, and without any heat source this gas would al- 
most instantly cool down to temperatures at which the Lya 
emissivity would be practically zero. Because the gas den- 
sities in the cold flows are high, this gas is capable of self- 
shielding and photoheating is negligible. We therefore re- 
sorted to 'gravitational heating', which refers to the hydro- 
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dynamical heating that the cold gas undergoes as it navi- 
gates down to the center of the dark matter potential. This 
gravitational heatin g mecha nism was described in an ana- 
lytic model by iHaiman et alj (i200Q 'l who predicted the ex- 
istence of Lya blo bs as a result of this cooling radiation. 
iFardal et al.l (|200ll 'l obtained similar results using hydrody- 
namical simulations -which already contained cold streams 
of gas, albeit at significantly lower resolution. In both pa- 
pers, the gas emitted all its gravitational binding energy as 
cooling radiation prior to reaching the central galaxy. This 
scenario would translate to /grav ^ 1 in our model. Our 
work may be viewed as an important improvement over - or 
an extension of - previous work, because our model is likely 
immune to various important objections that can be raised 
against older models: (i) we have shown that a significant 
of the gravitational binding energy may be dissipated in the 
hot gas without invalidating our conclusions (§ I4.6.2[) . {ii) 
in our model, gas accretion occurs along dense, compressed, 
flows that cover a tiny solid angle when viewed from the 
center of their host halo. Therefore, these cold flows could 
not be easily disrupted by powerful outflows from starburst 
activity or AGN, which would preferentially propagate into 
the lower density (volume filling) phase of the halo gas. This 
implies that cold flows can coexist with powerful sources , 
as observed for some of the LABs ([Chapman et al.l I2001I : 
iBasu-Zvch fc Schar j|2004l : iGeach et al.ll2007l , bOOgl TFor the 
same reason, we do expect these outflows to enrich the in- 
flowing cold streams with metals and dust, which could 
severely quench the Lya flux from the cold streams (see 

§113. 

More recently, iFurlanetto et al] (|2005h used high reso- 
lution SPH simulations to study Lya emission from struc- 
ture formation, which included cold flows. In their highest 
resolution simulations, IFurlanetto et al.l l|2005l ') have a spa- 
tial resolution of ~ 1 kpc which is comparable to that of 
iKeres et al.l l|2009l ). but in a cosmological volume that is too 
small to contain the massive halos that we associate with 
Lya blobs. Despite the absence of these massive halos, these 
simulations contained luminous blobs with properties that 
resemble the observed ones . Henc e, th ese models -just lik e 
the models of lHaiman et all (I2OO0I ') and lFardal et al] l|200ll )- 
overproduce the number density of blobs: existing observa- 
tions suggest that the field SSA21 contains a number density 
of blobs that is ~ 10— 10^ larger than the Universe as a whole 
(see § 13.21 and Fig HI . Lya blobs are likely rarer than was 
originally thought. 

Simulations thus seemed to create luminous Lya blobs 
(L-Lya ~ 10** erg s^^) around too abundant lower mass ha- 
los. This may partly be due to the fact that simulations 
contain a non-negligible fraction of cold gas that is locked 
up into denser clumps which reside both inside and outside 
the cold fiows. Because of their enhanced density, most of 
the Lya luminosity would come from these discrete clumps. 
These clumps may also be sites in which stars form which 
may boost the Lya emissivity from these clumped regions 
even more. However, it is likely that Lya emission does not 
escape efficiently from star forming regions ai z — 3: obser- 
vations of 2 = 3 Lyman break galaxies (LBGs) suggest that 
~ 20 — 25% of star forming galaxies has a strong enough 
observe d Lya emission lin e to classify as a Lya emitting 
galaxy l|Shaplev et al.|[2003l ). Furthermore, a significant frac- 
tion of these Lya emitting galaxies have emission lines that 



are weaker than expected on the basis of recombination the- 
ory (Dijkstra & Westra 2009). Without a good model for 
how Lya escapes from these clumps, their Lya luminos- 
ity may well have been overestimated. This could explain 
why simulations have overpredicted the number density of 
Lya emitting blobs. Al ternatively, the simula ted cold fiows 
are too warm (see e.g. IFurlanetto et al]|2005l . for a discus- 
sion on why gas temperatures in self-shielded regions are 
notoriously difficult to simulate). As was mentioned earlier 
I4.6.2|l . the gas temperature in our 'AMR model' was be- 
low that encountered in the simulation. This suggests that 
the AMR simulation also associates (much) more luminous 
Lya blobs to halos of a given halo mass than our model. 
This may result in a predicted number density of Lya blobs 
that is at odds with the observations. 



5 CONCLUSIONS 

Recent hydrodynamical simulations of the formation of 
galaxies show that baryons assemble into galaxies through 
a two-phase medium which contains filamentary streams 
of cold (T ~ 10* K ) gas in pressure equilibrium with a 
hot gaseous halo ( e.glKeres et a.l.ll200Ej: lOcvirk et al.ll2008l : 
iDekel et al.l l2009l : iKeres et all l2009l : Ugertz et al.l 120*091 ). 
These cold fiows c ontain ~ 5-25% of the total gas content 
(iKeres et all |2009|) i n halo s as massive as Mhaio ~ a few 
10"M^ (|Dekel et al.|[2009l ). 

At the typical densities and scales of the cold fiows 
{ric ^1 cm"''), the gas is self-shielded from external ionizing 
radiation (see §[3TT} a-nd is heated through gravitational con- 
traction. We have demonstrated that if >10% of the change 
in the gravitational binding energy of a cold flow goes into 
heating of the gas, then the simulated cold flows are spa- 
tially extended Lya sources with luminosities and number 
densities that are comparable to those of observed Lya blobs 
(see Fig|4)|. 

Furthermor e, the typical velocit y widths of the LABs 
span the range (|Matsuda et al.ll2006l ) of 500-1700 km s'^ , 
which is consistent with the model in which infall occurs 
at velocities in the range I.0-2.5ucirc ~ 220-1500 km s~^, 
which naturally result in the observed velocity widths (see 
ij |4.1l The filamentary structure of the cold flows may explain 
the wide range of observed Lya blob morphologies. 

The association with massive halos na t urally 
places LABs i n ov e rdens e regions (|Steidel et al. I I2OOOI : 
iMatsuda et al.l |2004 120061 ). Furthermore, the simulated 
cold flows are dense and cover only a small solid angle 
when viewed from the center of their host halo. Therefore, 
they could not be easily disrupted by powerful outflows 
from starburst activity or AGN, which would preferen- 
tially propagate into the lower density (volume filling) 
phase of the halo gas. This implies that cold flows can 
coexist with powerful sources, as observed for some of the 
LABs (IChapman et ahl I2OOII : iBasu-Zvch fc Scharj |2004| : 
iGeach et al.l |2007| , I2OO9I ). The association of LABs with 
sub-mm sources or AGN is not surprising, since these 
sources are triggered by the infall of cold gas into massive 
galaxies. Even when the central sources are sufficiently 
energetic to power the observed Lya emission, this causal 
relationship may be difficult to achieve (see § 14. 4p . Our 
model associates the spatially extended Lya emission 
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with the independent process of cooling within the dense 
streams of inflowing cold gas. This resolves the puzzling 
observation that very similar LABs ma y have completely 
different sources associated with them ( Geach et al.l 120071 : 
lYanget al.|[2009l ). 

Of course, it is possible that cold flows are partially 
photoionized by an associated h ard X-ray source (a s is ob- 
served for ~ 17% of the blobs 'Ge ach et af] 120091 '), or by 
star forming galaxies that are embedded within the cold 
flow (jFurlanetto et al.ll2005l l. In these cases, the heating rate 
would be (locally) boosted. This would only make our to- 
tal predicted Lya luminosities from the cold flows larger, 
and would strengthen our conclusion that Lya blobs are an 
observational signature of cold accretion into galaxies. 

Regardless of the precise heating mechanism, the cold 
flows should reveal their filamentary geometry in deeper and 
or higher resolution images (see 14. 5p . Once this gas geom- 
etry has been probed, it may be possible to constrain the 
dominant heating mechanisms of the cold gas as it navi- 
gates toward the center of the dark matter halo potential. 
Alternatively, if one does not find any evidence for filamen- 
tary spatially extended Lya sources in the near future, then 
this is equally interesting, because this may suggest that 
the simulated cold flows are different than the flows that 
occur in nature, or - worse - do not represent reality at all 
(see § 14.6.21 for a discussion on instabilities that may affect 
cold flow properties). In either case, Lya observations are 
expected to provide us with unique constraints on this in- 
triguing mode of gas accretion onto galaxies. 
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APPENDIX A: THE LYa BLOB LUMINOSITY 
FUNCTION IN AN OVERDENSE REGION 

ISteidel et alj (|l998l ) found an o verdensity of (5lbg = 5.0± 1.2 
Lyman Break Galaxies (LBGs [Steidel et al.ll2003l ) per unit 
volume within their survey volume of K ~ 0.15 x 10^ cMpc''. 
This highly overdense region is contained within the approx- 
imately nine times larger volume, Vs (the subscripts 's' and 
'S' refer to the smaller and larger survey v olumes, respec- 
tively), in which the Lya blobs were found (|Matsuda et al.l 
l2004il . In later work, iMatsuda et al.l (120051 ) found that the 
Lya blobs populate a region of our Universe that contains 
three fil aments. These filaments intersect right where Stei- 
del et al lSteidel et al.l (IT998 ) found their enhanced popula- 
tio n of LBGs. The r efore, the survey volume that was probed 
bv iMatsuda et all (|2004l ) is overdense, but no t as much as 
the 'su b-volum e' that was prob e d orig inally bv lSteidel et al.l 
(|2000l '). Indeed. IMatsuda et all l|2005l ) found that the num- 
ber density of Lya emitters in t he filaments - and therefo re 
the volume originally probed bv lSteidel et al.1 (|l998l . l2000l ')- 
is approximately two to three times higher than the average 
number density within their entire volume. If we assume that 
the LBG overdensity is suppressed by a similar fact or, then 
this im plies that the overall survey volume of .Matsuda et al.l 
(|2004l ') contains ~ (6.0± 1.2)/3 ~ 2.0 ±0.4 more LBGs than 
average. 

LBGs populate dark matte r halos of mass m — 
1O"'®M0 (| Adelberger et ahllioosl ). which are biased trac- 
ers of the overall mass density field. Therefore, we con- 
vert the LBG overdensity, 5lbg = 1-0 ± 0.4 into an over- 
all matter overdensity through <5lbg ~ b^BcSM, in which 
feLBG ~ 3.0 is the l inear bias parameter of dark matter halos 
hosting the LBGs (|Sheth et aLlf26oil '). We therefore find that 
5m = 0.3 ± 0.1, which corresponds to 5m = (3 ± 1) x cr(Vs). 
This result corresponds to the range of overdensities that 
was usecQ in this paper. 



^ iLidz et al.l ll2007t) compute the differential probability that a 
point is at an overdensity Sn when smoothed on a scale R, 
given that it is at an overdensity 5r when smoothed on a smaller 
scale r, P{5[i, aii\Sr, (^r) (see their Eq. 4). For a bias parameter 
''LBG = 3. we find that the m ass overdensity within the survey 
volume of ISteidel et al.l l ll998l) is 5m, r = = 1.7 ± 0.4 = 

(6.0 ± 1.5) X cr(V's). Here, the subscripts 'r' and 'R' refer to the 
smalle r and larger survey volumes, respectively. Using the expres- 
sion of lLidz et ahl ||2007|) we find that ^M.r = (6.0 ± 1.5) X a{Va) 
translates to i5m,r = (3.3±0.8) X a{Vs). This procedure therefore 
gives us an estimate of the overdensity that is consistent with our 
estimate based on the number of observed LBGs and Lya emit- 
ting galaxies. 



